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We report a measurement of the ortho-para transition 
rate in the p/j,p molecule. The experiment was conducted 
at TRIUMF via the measurement of the time dependence 
of the 5.2 MeV neutrons from muon capture in liquid hy- 
drogen. The measurement yielded an ortho-para rate A op = 
(11.1 ± 1.7±o:|) x 10 4 s -1 that is substantially larger than the 
earlier result of Bardin et al. We discuss the striking implica- 
tions for the proton's induced pseudoscalar coupling g p . 



The protons's weak axial current is governed by the 
weak axial form factor Fji(q ) and the induced pseu- 
doscalar form factor Fp(q 2 ) and their corresponding cou- 
pling constants g a — Fa(0) and g p = Fp(-0.88m 2 ). The 
values of the coupling constants and the q 2 dependence 
of the form factors are determined by the proton's sub- 
structure and the underlying QCD dynamics. The ap- 
proximate chiral symmetry of QCD implies a rigorous 
relationship between these constants that yields either 
g p = (6.5 ± 0.2) g a or g p = 8.2 ± 0.2 [1-4], and its ex- 
perimental verification is an important test of low energy 
QCD. 

The dominant uncertainty in extracting g p from /x 
capture in liquid H2 is the ortho-para transition rate 
A op from the ortho (i = 1) pfip state to the para 
(£ = 0) p/ip state. This transition has been inves- 
tigated theoretically by several authors [5-8]. It pro- 
ceeds by Auger emission and its rate is a function of 
the electronic environment of the pfip molecule, most no- 
tably analogues of the neutral Hi atom [(pfxp) + e\ and the 



charged H3 ion [(pfip) + i^2] + • The only published mea- 



surement of the rate (4.1 ± 1.4) x 10 



by Bardin et 



al. [9] is two standard deviations from the calculated rate 
(7.1 ± 1.1) x 10 4 s- 1 of Balakov et al. [6]. 

Given its importance in determining g p we have per- 
formed a new measurement of the rate A op . Both the 
Bardin et al. experiment [9] and the current experiment 
were based on the measurement of the time spectrum of 
the 5.2 MeV neutrons from /ip — > nv capture. Because 
of the spin dependence of the capture rate, and the spin 
configurations of the I = 0, 1 molecules, the rate A op is 
encoded in the neutron time dependence. Neglecting the 
initial formation of ortho molecules, the time dependence 



N(t) oc e — (Ap>jvf + (A OM - Ap^/)e op ), (1) 

where Ap> is the muon disappearance rate and Aom and 
Apm are the ortho and para capture rates. Since Ad 
and the hyperfine ratio Aom / Apm are well established 
by experiment or theory, the rate A op can be extracted 
from the neutron time spectrum. 

Our experiment was conducted on the M9B muon 
channel at the TRIUMF cyclotron. Incoming muons 
were counted in a two-element plastic scintillator beam 
telescope. The larger scintillator 51 (33x33x0.32 cm 3 ) 
registered particles emerging from the beam pipe and 
the smaller scintillator 52 (10x10x0.16 cm 3 ) registered 
particles incident on the H2 flask. The target [10,11] 
comprised a cylindrical flask of length 15 cm and di- 
ameter 16 cm containing 2.7 liters of liquid H2 with 
4±2 ppm deuterium contamination and <1 ppb contam- 
ination from Z>1 impurities. The target flask had Au 
front/side walls and a Cu back- plate and the surrounding 
vacuum vessel had cylindrical Ag walls and a mylar front 
window. 

Outgoing neutrons were detected in a liquid scintilla- 
tor counter array consisting of four counters (A2-A5) 
with diameters 12.7 cm and thicknesses 12.7 cm and one 
counter (Al) with a diameter 12.7 cm and a thickness 
5.1 cm. The scintillator materials were NE224 (Al), 
NE213 (A4) and BC501A (A2, A3, A5) which permit 
n/7 discrimination by pulse-shape. Sandwiched between 
the neutron counters and the vacuum vessel were plas- 
tic scintillator veto counters (E1-E5) for identification of 
/i-decay electrons. 

Signals from each neutron counter were fed to pulse- 
shape discrimination modules [12]. The modules gener- 
ated outputs indicating either a 7-ray trigger (G) or a 
neutron trigger (A) via the comparison of the integrated 
signals in a 420 ns prompt gate and a 52 ns delayed gate. 
To reduce the neutron trigger rate from mis-identified 7- 
rays/electrons we rejected those triggers accompanied by 
either a veto counter signal (E) or a large neutron counter 
signal (O). All events that fulfilled the neutron trigger 
(N-E-O) and pre-scaled events that fulfilled the 7/e trig- 
ger (G) were recorded. The pulse-shape data included 
both primary information from the pulse-shape modules 
(PSD1) and secondary information from a charge inte- 
grating ADC (PSD2). Histories of beam telescope hits 
and veto counter hits were recorded in a multihit TDC. 
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The experiment was conducted in four distinct data- 
sets and yielded a total of 3.7x 10 10 \i~ stops in liquid H 2 . 
Empty target runs (yielding 0.9 xlO 9 [i~ stops) and pos- 
itive muon runs (yielding 2.0 xlO 9 fi + stops) were used 
for background studies. In sorting the data we applied 
cuts to separate the neutrons from /up capture, with en- 
ergies 5.2 MeV and lifetimes ^2 /us, from backgrounds 
which included both mis-identified 7-rays from electron 
bremsstrahlung and neutrons from nuclear /1 capture, 
photo-nuclear reactions, accelerator sources, cosmic-rays, 
etc. 

A pulse-shape cut and a pulse pile-up cut were used to 
minimize the contamination from mis-identified 7-rays. 
The pulse-shape cut involved a comparison of the inte- 
grated amplitudes in the signal peak and the signal tail, 
and used both PSD1 and PSD2 to achieve optimal n/7 
discrimination. The pulse pile-up cut rejected any event 
with >1 neutron counter hits in a preceding 600 ns time 
window, such pile-up potentially distorting the n/7 dis- 
crimination. 

An energy cut was used to enhance the signal of recoil 
protons from 5.2 MeV neutrons by requiring an energy 
equivalent of typically 2.5-5.0 MeV. A blank cut was used 
to reject the events from fi capture in high-Z materials 
(i.e. tau = 73 ns, ta 9 = 87 ns and Tc u = 163 ns [13]) by 
rejecting any event with >1 SI hits in a 400 ns window 
preceding the neutron. Finally, we used a duty cycle cut 
to reject any events occurring within 32 fis of the beam- 
off period in the cyclotron time structure, and an errant 
muon cut to reject any events occurring within 5 /is of 
an incoming beam particle that fired SI but missed 52. 

The time spectra between incoming muons and de- 
tected neutrons for events that survived all cuts exclud- 
ing the blank cut and all cuts including the blank cut 
are plotted in Fig. 1. Note that for each event we in- 
crement the time spectrum for every S1-S2 coincidence 
in the multihit TDC. The spectra show a short-lifetime 
(r ~ 80 ns) component that is most evident for events 
that failed the blank cut, a long- lifetime (r ~ 2 /is) com- 
ponent that is only evident for events that passed the 
blank cut, and a continuum background. The short- 
lifetime component is dominated by Ag/Au capture, the 
long-lifetime component is dominated by H2 capture, and 
the continuum arises from uncorrelated muons. 

The procedures for accounting for backgrounds in the 
time spectrum of the events passing all cuts were (i) sub- 
traction of the continuum background from uncorrelated 
muons, (ii) determination of contributions from photo- 
nuclear reactions following electron bremsstrahlung, and 
(iii) the determination of contributions from /1 capture in 
low-Z materials. 

The continuum is not perfectly time independent as de- 
cay electrons from uncorrelated muons can reject events 
by firing the veto counter (E) or overload circuit (O). 
This random vetoing of neutron events by uncorrelated 
muons can be seen as the shallow hole for t < 4 fis 



for events failing the blank cut. Note that for times 
t >> TAu/Ag the events failing the blank cut are en- 
tirely dominated by the continuum background while the 
events passing the blank cut are comprised of /up sig- 
nal and continuum background. Thus, by normalizing 
and subtracting the events that failed the blank cut from 
those that passed the blank cut, the continuum is re- 
moved. 
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FIG. 1. Time spectra between incoming muons and de- 
tected neutrons for events that passed all cuts except the blank 
cut (dashed line) and all cuts including the blank cut (solid 
line). The shallow hole for t < 4 us that's visible for events 
failing the blank cut is due to random vetoing. Note that 
t = —30 to 370 ns corresponds to the blank cut and the 
dashed and solid curves are normalized at t > 10 fis region. 

The photo-nuclear background originates from the 
photo-nuclear interactions [e.g. 12 C(7,p)] in the liq- 
uid scintillators. Consequently, a tiny fraction of 
bremsstrahlung photons from decay electrons yield pro- 
ton recoils. This background carries the 2.2 fj,s lifetime of 
muonic H2 and therefore distorts the effective disappear- 
ance rate of /ip capture neutrons. To measure this back- 
ground we conducted a fi + measurement where neutrons 
from [i capture are absent but those from photo-nuclear 
reactions are present. From the amplitude of the 2.2 [is 
lifetime component in the \x + neutron time spectrum we 
determined a contribution of photo-nuclear interactions 
to /i~ running that averaged about 14% for the thick 
counters N2-N5 but rose to 55±8% for the thin counter 
ATI. We used a GEANT3 simulation [14] to account for 
the difference in the 7-ray spectra from /i + and \i~ stops 
due to e + annihilation. 

The background from ji capture in low-Z materials, e.g. 
scintillators, light-pipes, is worrisome as the muonic C 
and H2 lifetimes are similar (2.028 /us and 2.195 fjs). This 
background was determined by performing an empty tar- 
get run in which fjT stops occur in the Cu back-plate not 
the liquid H2. A GEANT3 simulation was used to es- 
tablish that stops in scintillators, light guides, etc., were 
dominated by scattering from the beam pipe window and 
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the beam telescope and thus unaffected by the target 
filling. The measurement showed no evidence for back- 
grounds from low-Z capture and established a limit of 
<7% on carbon backgrounds in production running. 
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FIG. 2. Time spectrum of events that survive all cuts after 
the continuum subtraction. The points are the data and the 
solid curve is the benchmark fit. 

The time spectrum of neutron events that survived all 
cuts, after the subtraction of the continuum, is plotted 
in Fig. 2. A typical energy spectrum (JV4, data-set I) 
is plotted in Fig. 3. It shows the expected spectrum for 
5.2 MeV neutrons as obtained by a GEANT3 simulation 
of the experimental set-up. 

When muons stop in liquid H2 they evolve through a 
sequence of different atomic and molecular states. Also, 
if D 2 is present other species of /i-atoms and /i-molecules 
are generated [15,16]. Therefore corrections were applied 
to Eqn. 1 to account for the p/ip formation from the sin- 
glet /ip atoms and the p/id formation due to D2 contam- 
ination. We assumed that, due to their typically lower 
energies, the efficiency for detecting neutrons from d/ 3 He 
capture is significantly smaller than the efficiency for de- 
tecting neutrons from proton capture, this making proton 
capture in p/id molecules the dominant correction from 
D2 contamination. 

The determination of A op was obtained by fitting to 
the Fig. 2 time spectrum. The fitting function was based 
on Eqn. 1, but supplemented with corrections account- 
ing for p/ip formation and D2 contamination. Wc also 
included a correction term with disappearance rate Ad 
to account for photo-nuclear backgrounds. In most fits 
we fixed the ratio Aqu/^-PM — 2.59 [2,3] and rate 
Ad = 4.55 x 10 5 s -1 [19] at their known values and only 
varied the overall normalization and rate A op . 

Our benchmark fit (fit A) to all counters and all data- 
sets, which demonstrates the good agreement between 
the measured spectrum and the theoretical spectrum, is 
shown in Fig. 2. This yielded a rate A op = (11.1 ± 1.7) x 
10 4 s~ 4 with a chi-squared p.d.f. of 1.03 for a fit range 



of t = 1.3 to 24.3 /is, i.e. beginning beyond the region 
contaminated by Au/Ag capture. 

The results of fits to study the sensitivity to back- 
grounds and corrections are listed in Table I. It shows 
that omitting the p/ip formation correction (fit B) or 
p/id contamination (fit C) correction has little effect 
(< 1 a) on A op . Including a 7% carbon background 
(fit D) changed the rate from (11.1 ± 1.7) x 10 4 s" 1 to 
(11.7 ± 2.0) x 10 4 s _1 , a 0.5 a effect, while omitting the 
photo-nuclear background (fit E) changed the rate from 
(11.1 ± 1.7) x 10 4 s" 1 to (8.0 ± 1.1) x 10 4 s~\ a 1.8 a 
effect. Other tests indicated that varying the hyperfine 
ratio Aqm I Apu by ±5% changed the ortho-para rate 
A op by ±4%. 
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FIG. 3. Energy spectrum for counter N4 and data-set I 
for events surviving all cuts with times t = 1.1-3.3 fis after 
continuum subtraction. The points are the data and the his- 
togram is the simulation The data's upper edge (chan. ~550) 
corresponds to 5.2 MeV proton recoils and the data's lower 
edge (chan. ~250) results from the pulse-shape cut. 

Also given in Table I are results from fits to the in- 
dividual time spectra of the five neutron counters, thus 
demonstrating the consistency between counters. Similar 
checks of separate fits to the four different data-sets from 
the two running periods, which involved some differences 
in calibrations and settings, yielded results ranging from 
A op = (10.3 ± 2.6) x 10 4 s- 1 to (13.7 ± 3.0) x 10 4 s" 1 . 
Lastly, we found no evidence of significant sensitivity to 
either the start-time or the end-time of the fit. 

Several other tests of data integrity were performed. 
We sub-divided the neutron events into a low energy 
data-set (typically 2.5-3.7 MeV) and a high energy data- 
set (typically 3.7-5.0 MeV) and found consistent results 
for A op . We analyzed the neutron events with energies 
> 5.2 MeV and found no evidence for a long- lifetime 
component. Finally, the time spectrum of electron events 
yielded r = 2.198 ± 0.002 /is and the time spectrum of 
empty-target neutron events yielded r = 163 ± 1 ns, con- 
sistent with the known values for the /ip lifetime [19] and 
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the nCu lifetime [13]. 

Our final result is A op = (11.1 ± 1.7±g;§) x 10 4 s" 1 
where the first error is statistical and the second error 
accounts for uncertainties in background determinations, 
chemistry corrections and Aom / Apm (±5%). Our result 
is somewhat larger (by 1.9 a) than the theoretical esti- 
mate A op = (7.1 ± 1.1) x 10 5 s" 1 of Bakalov et al. [6] 
and substantially larger (by 3.1 a) than the earlier mea- 
surement A op = (4.1 ± 1.4) x 10 4 s _1 of Bardin et al. 
[9]. 

In extracting g p from /i capture in liquid H2 the differ- 
ence between A op = 4.1 x 10 4 s _1 and 11.1 x 10 4 s _1 is 
striking. With the Bardin et al. value A op = 4.1 x 10 4 s" 1 
the TRIUMF RMC experiment [17,18] yields g p = 12. 2± 
1.1 [17,18], the Saclay OMC experiment [19] yields g p = 
10.6 ± 2.7 [19], and the average value of the earlier liq- 
uid H 2 experiments is g p = 10.4 ±4.1 [20,21], with the 
TRIUMF result being clearly inconsistent with theory. 
With the new measured value A op = 11.1 x 10 4 s _1 the 
TRIUMF RMC experiment yields g p = 10.6 ± 1.1, the 
Saclay OMC experiment yields g p — 0.8 ± 2.7, and the 
average value of the earlier liquid H2 experiments yields 
g p = 5.6 ±4.1, with the Saclay result being clearly incon- 
sistent with theory. 

The inconsistency between the present value for A op 
and that measured by Bardin et al. is disconcerting. Note 
that an important difference between the two experi- 
ments is the beam time structure, the earlier work using 
a pulsed beam while the current work using a continuous 
beam. The experiments also differed in the total num- 
ber of 5.2 MeV neutrons, the signal-to-noise for 5.2 MeV 
neutrons, and the presence in the Bardin et al. experi- 
ment of a long-lifetime background. However, we see no 
obvious explanation for the conflicting results. 

Note that one consequence of the different time struc- 
tures in the two experiments is the typical time between 
the parent muon and the detected neutron, in the rele- 
vant region of the time spectra, is greater in the Bardin et 
al. experiment than the current experiment. Therefore, 
at least in principle, it is possible that unconventional 
/i-chemistry in liquid hydrogen could manifest itself as 
different neutron disappearance rates in the two experi- 
ments. 

In summary, we report a determination of the pp,p 
ortho-para rate from the time spectrum of the 5.2 MeV 
HP capture neutrons. Our result A op = (11.1 ± 1.7±o'§) x 
10 4 s _1 is somewhat larger than the theoretical result of 
Bakalov et al. [6] and strikingly larger than the only pre- 
viously published measurement [9]. The new result has 
dramatic consequences for the determination of the in- 
duced pseudoscalar coupling g p . 
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TABLE I. Results for the rate A op from the fits to the con- 
tinuum-subtracted time spectra. Values are given for both the 
summed spectrum of all neutron counters and the individual 
spectra of each neutron counter. See text for explanation of 
fits A,B, C, D and E. The units of A op are xlO 4 s -1 . 

Cntr. Fit A. Fit B Fit C Fit D Fit E 

~AU 11.1±1.7 12.2±1.8 11.6±1.7 11.7±2.0 8.0±1.1 

Nl 7.2±10.6 8.2±11.1 7.7±10.6 7.3±10.9 2.2±3.9 

N2 13.8±5.5 15.4±5.6 14.4±5.5 14.8±6.5 8.3±3.1 

N3 8.9±2.1 9.9±2.3 9.4±2.2 9.2±2.3 7.1±1.7 

N4 14.6±6.2 16.1±6.7 15.1±6.2 15.5±7.3 10.3±4.1 

N5 12.2±3.0 13.5±3.0 12.7±3.0 12.9±3.4 10.6±2.6 
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